Saddle point states in 2D superconducting film biased near the depairing current 



(N 

O 

(N 



D.Y. Vodolazo\0 

Institute for Physics of Microstructures, Russian Academy of Sciences, 603950, Nizhny Novgorod, GSP-105, Russia 

(Dated; March 2, 2013) 

The structure and energy of saddle point (SP) states in 2D superconducting film of finite width 
w with transport current I are found in the framework of Ginzburg-Landau model. We show that 
very near depairing current Idep the SP state with a vortex does not exist and it transforms to 2D 
nucleus state, which is a finite region with partially suppressed order parameter. It is also shown 
that for slightly lower currents the contribution of the vortex core energy is important for SP state 
with a vortex and it cannot be neglected for I > 0.6/dep. It is demonstrated that in the film with 
local current concentration near the bend the energy of SP state may be much smaller than one in 
the straight film and it favors the effect of fluctuations in such a samples. 

PACS numbers: 74.25.0p, 74.20.De, 73.23.-b 
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I. INTRODUCTION 

Many physical systems have several metastable states 
at fixed external parameters (temperature, magnetic 
field, etc.) which correspond to different local minima of 
their free energy. As an example it could be mentioned 
the different configuration of DNA molecule [l|, vortex 
'molecules' in the mesoscopic superconductors 0, [1], or 
magnetic states of the lattice of magnetic nanocaps [I]- 
The metastable states are usually separated from each 
other by the energy barrier AF which could be overcome 
due to thermal activation. If the height of the energy 
barrier is much larger than the thermal energy ksT then 
such a transitions are relatively rare events and one may 
estimate their rate as ^ exp{—AF/kBT) Q. The stan- 
dard procedure to calculate AF is to find the saddle point 
states which correspond to the local maxima of the free 
energy and find the energy difference between them and 
the metastable states [ll-ll, Q . 

In the superconductors with transport current less 
than critical one this idealogy could be also applied for 
calculation of the finite resistance appearing due to fluc- 
tuations. For ID superconducting wires Langer and Am- 
begaokar (LA) Q in framework of the Ginzburg-Landau 
(GL) model find that in the saddle point state the su- 
perconducting order parameter ^p = {iple^^ is partially 
suppressed in the finite region along a wire and the am- 
plitude of suppression depends on the current. If one 
starts from such a ID nucleus state the time evolution 
of the order parameter will inevitably lead to the phase 
slip [6,] and to the finite resistance Rwire- LA find depen- 
dence of AF both on the current and temperature which 
with good accuracy Q could be written as 
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where Fq ^ ^^d / Wtt^ , $0 is a magnetic flux quantum, 



w is a width and d is a thickness of the wire respec- 
tively, A is a London penetration depth, ^ is a coherence 
length and Idep = c^oU!d/12\/3n^^X^ is a depairing cur- 
rent in the Ginzburg-Landau model. According to the 
general concept Rwire = v exp{—AF]^A/kBT) where the 
pre-exponential factor v is calculated in Refs. tS, i^- 

In contrast to ID wire in thin [d ^ A) 2D film several 
saddle point states exist at given value of the current. 
Further we discuss the case of the relatively narrow film 
with £^ w < X^/d (this condition ensures the unifor- 
mity of the current distribution over the width of the su- 
perconductor with transport current in the ground state). 
At the present time it is distinguished three types of the 
saddle point states in such a samples: i) the state with a 
de vortex (Tol - [l5| . ii) the vortex- antivortex (VA) s tate 
■[17 [and iii) LA state extended to 2D case |13l - [l5l |. In 
Ref. [l5j| it was argued that VA state has at least twice 
larger energy than the single vortex state and LA-like 
state is the most energetically unfavorable among the 
considered states at all currents / < Idep- However in 
Ref. [15| the current dependence of AFla (see Eq. (1)) 
was not taken into account. If one takes the result for 
the energy of the single vortex (V) state near depairing 
current 1 — 1/ Idep ^ 1 (found in the London model [Tsl - 
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and compare it with Eq. (1) then it is easy to see that 
even for wide films w 3> ^ there is a finite region of the 
currents 1 - {3tt^/2V2w)'^ < I /Idep < 1 where AFla < 
AFy. But such a quantitative comparison is not correct 
because Eq. (2) does not contain the energy of the vortex 
core {Ecore ~ 0.38Fo for vortex located far from the edges 
p^ ) and, as we show below, it seriously underestimates 

AFy at / > 0.6/dep- 

The effect of the vortex core was taken into account 
in the recent work [l^] using GL model. But authors fo- 
cused on the region of small currents / ^ Idep (where 
contribution of Ecore to AFy is relatively small) and 
they concluded that for films with w > Wc — 4.4^ the 
energy of single vortex state is lower than the energy of 
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LA state. Below we show that vortex SP state does not 
exist at / ^ Idep and it transforms to vortex free 2D nu- 
cleus (which is a finite region with partially suppressed 
order parameter) located near the edge of the film. We 
confirm our above estimation that very near Idep even 
for wide films w ^ Wc saddle point state of Langer and 
Ambegaokar has a lowest energy among over SP states. 
We also study practically important case of the film with 
one 180° bend which models the part of the supercon- 
ducting meander used in superconducting single photon 
detectors [l^-[2l|. We show that due to current concen- 
tration near the bend the jump to the saddle point state 
needs much less energy at J — >■ /c in comparison with 2D 
film with uniform current distribution, which promotes 
the effect of fiuctuations in such a samples. This result 
is rather general and could be applied to any supercon- 
ducting system with a local current concentration. 

II. METHOD 

To find the state which is a solution of stationary 
Ginzburg-Landau equations but is unstable one (saddle 
point state) we use the following numerical procedure. 
In case when we search for SP state with a vortex we 
first put (as an initial condition) the phase distribution 
corresponding to the vortex seating in the point (n,m) of 
the discrete grid and additionally fix the phase difference 
between adjacent points ip{n, m + 1) — fin, m — 1) = tt 
(which, in some respect, pins the vortex at the point 
(n,m)) at any time step. The solution of the stationary 
GL equations we find by using the relaxation method 
(by adding the time derivative dip/dt in GL equations 
and waiting when it goes to zero). By variation of the 
current (as an external parameter) it is possible to find 
a stationary state when the vortex does not move from 
the point (n,m). At low currents there are several points 
where the vortex position could be fixed by this method 
for given /. By our definition the one with the highest 
energy corresponds to SP state (it corresponds to the lo- 
cal maximum of Gibbs energy as function of the vortex 
position in the London approach [iS-E [3 ) ■ After 
reaching such a state the energy difference could be found 
using the following expression 

Ai^ = Fsaddle — Fground — 2^^^^ ('^) 

where Atp is an additional phase difference between ends 
of the film which appears in SP state in comparison 
with the ground state and Fsaddle and Fground are the 
Ginzburg-Landau free energies of the saddle point and 
ground states respectively. 

To find the vortex free saddle point state we fix the 
magnitude of the order parameter \^\ > in one point 
at the edge and allow to vary -0 in all other points. Than 
we increase the current up to the moment when such a 
state becomes nonstationary. By our definition we find 
the vortex free SP state corresponding to the given value 



of the current. We checked that if we start from this 
state and let if) vary everywhere in the film, the vortex 
is nucleated in the point where we initially fix |^| and 
passes across the film. This finding is an extension to 2D 
case of the main idea of LA that if one starts from SP 
state with finite everywhere in the sample the time 
evolution of the order parameter will inevitably lead to 
the phase slip and voltage pulse. 

The proposed method is much simpler (from point of 
view numerical procedure) than the methods used in 
Refs. [12, [2^ for finding SP states in the mesoscopic 
superconducting disk in magnetic field or 2D film with 
a transport current [ij, Similar procedure to fix 

the vortex position was used, for example, in Ref. [23 |. 
Moreover our method could be easily applied to 2D sam- 
ples of arbitrary geometry (triangles, disks, stars, etc.) 
and to 3D case to find the vortex free saddle point states. 
We checked that the method gives the same results for 
ID saddle point state (both for the spatial dependence of 
order parameter and excess energy AF) as analytically 
found in Ref. [1]. Validity of the proposed method is 
also supported by our results for 2D film because they 
coincide with results found in the London model at low 
currents. 

In numerical calculations the step of the rectangular 
grid was 5x = 5y = S^/A and width of the film varied 
from 4.5 up to 30 The length of the film was chosen 
L = Aw (which is long enough to neglect the effect of 
finite length). The boundary conditions, written here in 
dimensionless units (for units see for example [23!|), were 
the following: 7/'*V7/'|y^±/^/2 = il/wd and yip\x=±w/2 = 
0. 



III. RESULTS 

In Fig. 1 the dependence of AF on the current for 
different SP states is present for the film with w — 4.5^. 
We should note that vortex in such a narrow film has 
strongly deformed core (see upper inset in Fig. 1) and 
it resembles more a Josephson vortex than an Abrikosov 
one [2^. In wider films (see Figs. 2,3) the deformation 
of a vortex core occurs when the vortex seats near the 
edge on the distance Ax < 2^ (see insets in Figs. 2,3). 
Similar result was found for the vortex placed near the 
artificial defect (see Fig. 2 in Ref. [l^l) and the edge of 
the superconducting disk (see Fig. 3 in Ref. or film 
(see Fig. 2 in Ref. [l3). 

Unfortunately our numerical method does not allow to 
find SP state with a vortex when Ax < 1.5^ (last red 
circle at the largest current in Figs. (1-3) corresponds 
to Ax = 1.5^) because we were not able to find the sta- 
tionary solution of Ginzburg-Landau equation with ad- 
ditional condition ip{n,m -I- 1) — ip(n,m — 1) ~ t: (this 
effective pinning 'force' becomes not sufficiently strong 
to pin the vortex very near the edge). But we notice 
that if we fix \tp\ = along finite length line near the 
edge (see inset in Fig. 1 where contour plot of \ip\ at 




FIG. 1: Energy of saddle point states of three kinds: LA 
(green dashed curve), vortex/'line' (red circles/rare empty 
squares) and 2D nucleus (dense empty squares) in the film 
with w = 4.5^. In the insets we present the contour plot of 
at different currents and different kinds of SP states. 



I / Idep — 0.51 is shown) and find the stationary solution 
of GL equations with this additional condition then the 
excess energy AF of such a 'line' state (see rare empty 
squares in Figs. 1-3) is close to the energy of the vortex 
state (when vortex seats relatively close to the edge - see 
Figs. 2,3). Because of that one may approximate the 
energy of the vortex state when Ax < 1.5^ by the energy 
of the 'line' state with length I < 2£^. 

At currents / ~ Idep the vortex('line') state transforms 
to the vortex-free SP state (dense empty squares in Figs. 
1-3) when the phase circulation along any closed contour 
in the film is equal to zero and |'0| > everywhere in the 
film. To find it we fix the amplitude of the order param- 
eter j-^l in one point at the edge. Because of proximity 
effect and that / ~ Idep it leads to suppression oi \tp\ in 
the relatively large region around this point (see inset in 
Fig. 1 at I /Idep = 0.94). Further we call it as 2D nu- 
cleus state to distinguish it from ID nucleus state of LA 
(compare insets in Fig. 1 at / = and at I /Idep = 0.94). 

In Figs. (2,3) we also plot dependence of AF for the 
energy of the vortex SP state found in the London limit 
[l3l - [l5| (solid blue curve) 




where a = 3-\/37r^/4w and we add the energy of the vor- 
tex core Eeore = £^0- Numerical coefficient e is found 
from comparison of Eq. (4) with numerical result at 
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FIG. 2: Energy of saddle point states of three kinds: LA 
(green dashed curve), vortex/'line' (red circles/rare empty 
squares) and 2D nucleus (dense empty squares) in the film 
with w = 7^. Blue curve corresponds to Eq. (4). Black dot- 
ted line is the dependence AF/Fq = 1.43(1 - //1.026/dep) 
which fits well (deviation less than 2%) our numerical results 
at 0.6 < I/Idep < 0.97. At 7 > 7* ~ 0.737^^^ LA state has a 
lowest energy. 



TABLE I: Values of coefficients in the fitting expressions for 
energy of vortex/2D nucleus saddle point states (Eq. (5)) and 
vortex core energy when vortex seats in the center of the film 
(coefficient e in Eq. (4)). 



w/^ 


A 


B 


C 


n 


e 


7 


1.43 


1.026 


0.89 


0.7 


0.37 


10 


1.67 


1.026 


1.02 


0.7 


0.38 


15 


1.85 


1.028 


0.88 


0.6 


0.38 


30 


1.88 


1.034 


0.68 


0.5 


0.38 



7 = and it is present in Table 1 for different widths. No- 
tice the good agreement between GL and London model 
at / < O.dldep- At larger currents vortex is located on 
the distance Ax < 2^ from the edge and one has to take 
into account deformation of the core which provides the 
dependence e(7). It leads to large discrepancy between 
London (with e — const) and GL models at 7 > O.GIdep 
(see Figs. 2,3). Moreover at 7 ^ Idep vortex state trans- 
forms to the 2D nucleus state which cannot be found in 
the London limit. 

Our numerical results at I /Idep ^ 0.6 could be fitted 
(examples of fitting are present in Fig. 2,3 and inset in 
Fig. 5) by the following functions 

AF ( A{1- I/BIdep), 0.6 < I/Idep < 0.97, , . 

Fo ~\ C{1 - I/IdepT, 0.95 < I/Idep < 1. ^ ' 

Coefficients A,B,C and power n for different widths are 
listed in Table 1. Note that coefficient A is almost twice 



4 



larger than the result which follows from the London 
model (see Eq. (2)). It reflects the contribution of 
Ecore{I) to AF which for / > 0.6/dep cannot be neglected 
and Ecore gradually decreases with increasing /. 
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vortex state (London model)_ 
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FIG. 3: Energy of saddle point states of three kinds: LA 
(green dashed curve), vortex/'line' (red circles/rare empty 
squares) and 2D nucleus (dense empty squares near depair- 
ing current) in the film with w = 15^. Blue curve corre- 
sponds to Eq. (4). Black dotted line is the dependence 
AF/Fo = 1.85(1 - //1.028/dep) which fits well (deviation less 
than 2%) our numerical results at 0.6 < I/Idep < 0.97. At 
/ > /* ~ 0.92Idep LA state has a lowest energy. 

At I/Idep ^ 1 the power n < 1 (see Table 1) in Eq. (5) 
tells one that there is a finite (but rather narrow for wide 
films with w ^ ^) interval of currents I*{w) < I < Idep 
where ID nucleus (LA) state has the lowest energy (see 
Figs. 2,3). This counterintuitive, at first sight, result 
is explained by the presence of the last term in right 
hand side of Eq. (3). Although in LA state the order 
parameter is suppressed over whole width (see inset at 
/ = in Fig. 1) and it costs more condensation energy 
than in vortex/2D nucleus state, the phase difference Aip 
is much larger in LA state than in other SP states at 
I ^ Idep and it brings above result. 

Before we consider only single vortex state and 2D nu- 
cleus which is located near the edge of the film. In Fig. 
4 we demonstrate that the energy of 2D nucleus located 
in the center of the film is larger than the energy of the 
edge 2D nucleus. The vortex/antivortex state which is 
transformed from 2D nucleus state at lower currents (see 
solid and rare empty squares in Fig. 4) has energy larger 
than the single vortex state. The similar result was found 
in London limit where the difference reaches two times 
between AFy and AFva [11113 ■ 

We also find the saddle point states in 2D supercon- 
ducting film with a 180° bend - see left inset in Fig. 5. In 
our simulations we choose the width of the film w = 10^, 




w=15^ 
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FIG. 4: Energy of saddle point states of different kinds: sin- 
gle vortex (red circles), vortex-antivortex pair (blue squares) 
and 'line'/2D nucleus located at the edge (empty rare/dense 
circles) and in the center (empty rare/dense squares) of the 
film. 



the length of the sample in the bend region L=2w and 
two widths of the slit: Wsut — 2.5^ (shown in the left 
inset in Fig. 5) and Wsiu ~ 10^. In Fig. 5 we present our 
results for the energy of single vortex and edge 2D nu- 
cleus states. Note that the current in Fig. 5 is normalized 
to the critical current of the sample and not to the Idep 
as in Figs. 1-4 (Jc = 0.85Idep for film with Wsut — 2.5^, 



0.91 Idep for film with w. 



lit 



10^ and Ic 



^dep 



for 



straight film). 

We want to stress here that at I ^ Ic the energy of SP 
states is considerably lower in the film with bend than 
without it (compare results present in right inset in Fig. 
5) taken at the same ratio I / Ic- And the stronger the cur- 
rent concentration near the bend (which is manifested in 
lower value of Ic for smaller Wsnt) the smaller AF. We 
expect that the smallest AF could be reached in case of 
infinitesimally narrow and long crack near the edge of the 
film (26l . [2^ which provides the maximal current concen- 
tration and maximal suppression of Ic- We explain this 
effect by the partial suppression of the order parameter 
(on the scale of about several ^) in the region with the 
strongest current concentration even in the ground state. 
As a result it takes less energy to jump to SP state from 
the ground state with already locally suppressed {tpl- The 
proof of this idea also comes from our results at low cur- 
rents where AF differs a little for straight film and film 
with a bend (see Fig. 5). At low currents the suppres- 
sion of |^/;| near the bend is weak and it slightly infiuences 
AF. 

The decrease of AF in the ID wire with constriction 
was also found in Ref. [2^. But authors of [l^ did not 
present AF as function oi I /Ic and it is not clear is their 
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result connected with reduction of Ic (which gives trivial 
AF — ^ at current I = Ic < Idep and AF < AF^a at 
/ < Ic) or with a change in the power n in asymptotics 
AF ^ (1 — ///c)" as in our case (see right inset in Fig. 
5). 




0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 5: Energy of saddle point states (vortex and 2D nucleus) 
for the film with a bend and slits of different width: Wsut = 
2.5^ (red circles) and Wsut = 10^ (blue triangles) and the 
film without bend with w — 10^ (black squares). In the right 
inset we present the zoom at J ~ Jc which shows the energy 
of 2D nucleus SP states with dashed curves as the fitting of 
the numerical results. Green dashed curve shows the energy 
of LA state for straight film with w = 10^. 



for wide films w ^ £,) the lowest energy SP state corre- 
sponds either to edge 2D nucleus state or to state with 
a vortex located next to the edge and having strongly 
modified core. At currents / < O.GIdep (for films with 
w > 7^) the lowest energy saddle point state corresponds 
to the single vortex with ordinary core and the results of 
the London model are recovered. We demonstrate that 
in the film with current concentration, which leads to 
the local spatial variation of the order parameter in the 
ground state, the energy of 2D nucleus SP state may be 
much lower than the energy of any saddle point states in 
the film with uniform current distribution taken at the 
same ratio I / Ic ^ 1- 

The last result has an important practical consequence. 
It shows that if in the sample there are places with strong 
current concentration (bends, geometrical defects of the 
edge and so on) it favors the effect of fluctuations near 
Ic because of much low, in comparison with straight film 
without defects, energy barrier. For example, if ksT = 
O.lFo the switching of the straight film with w — 10^ to 
resistive state roughly occurs at / ~ 0.96/c = 0.96Idep 
(at this current AF = O.IFq - see right inset in Fig. 5). 
But for the film with a bend and Wsut = 2.5^ the energy 
barrier AF = O.lFo is reached at / ~ 0.85/c ~ 0.72Idep 
and transition to resistive state occurs at much lower 
current than one could expect from results for straight 
film. 

We believe that the found results are valid not only 
near Tc (where Ginzburg-Landau model is quantitatively 
valid) but at lower temperatures too. Indeed, in the re- 
cent work [29] asymptotics AF ^ (1 — I/Idep)^^'^ for 
/ — > Idcp was confirmed for LA like state in 'dirty' ID 
superconducting wire at any temperature. In the 'pure' 
ID wire only 15% difference with LA result at T — > 
was noticed in Ref. [SOj] in the limit / 31]. 



IV. CONCLUSION 

It is found that in 2D superconducting film with spa- 
tially uniform current distribution in the ground state 
the lowest energy saddle point state at J Idep corre- 
sponds to the state of Langer and Ambegaokar found 
for ID wire. Only at / < /* (where I*{'w) < Idep even 
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